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The recent real time KARMEN anomaly for the electron neutrino counting rates (obtained from stopped 
muon decays) is analyzed by employing a neutrino flavor rotation model in which the muon neutrino is in a 
superposition of only two mass eigenstates. On the basis of experimental electron neutrino counting oscillations, 
we find a neutrino mass splitting \m\ - rn{\ = (0.22 ± 0.02) (eU/c 2 ) 2 and a flavor rotation angle <j> = 0.34 ± 0.10, 
both within a 95% confidence interval. 
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Recent real time observations of the electron 
neutrino i> e , produced from a /i + decay 

M + -> e+ + + (1) 

indicate an oscillatory counting ratc[|lj. These 
data provide experimental evidence for a neutrino 
mass matrix with flavor rotations. The oscilla- 
tions of interest are made manifest in the experi- 
mental KARMEN anomaly 0. Our purpose is to 
discuss this fact, employing a model in which the 
muon neutrino is in a superposition of two mass 
eigenstates 

<v v \ = cos <j) < v x | + sin <f> < v-i \ , (2) 

with mass eigenvalues of mj and 7712, respectively. 

One may view Eq.(l) as occurring in two 
stages: (i) 

M+ -> W+ ff + (3) 
and (ii) 

W+ f ^e + +u e . (4) 

The virtual "effective" W%f standard elec- 
troweak charged Boson is far off the M-w mass 
shell, and appears in the lowest order tree dia- 
gram for Eq.(l) via the charged Boson propagator 
T> a p. The effective wave function is given by 

Wt ff ,M = J V aP (x - y)J + ' /3 (y)d i y, (5) 



where the charged current J + is the source of 
the W^fj Boson. In the far off M\y mass shell 
regime of Eq.(l), it is sufficiently accurate to write 
the wave function in Eq.(5) using the field-current 
identity of the Fermi formulation 

w <" = 9w ^S^ a{i ~ i5)N ^ (6) 

Eq.(6) employs units for which the W fine struc- 
ture constant aw = (9w/^ c ) ^ s related to 
the Fermi coupling fine structure constant by 
(V2)a w = (G F M W /Hc). 

In Eq.(6), the muon wave function ^ is the one 
particle to vacuum matrix element of the muon 
field operator ip. The muon wave function is 

*(x) =< 0\ip(x)\n + > . (7) 

The muon neutrino wave function iV M is the vac- 
uum to one particle matrix element of the muon 
neutrino field operator, 

Vfj,(x) — cos4>ui(x) + sintfiv2(x). (8) 

The neutrino wave function is 

N„{x) =< fyMaOlO > . (9) 

From Eqs.(2), (8) and (9), one finds that the 
muon neutrino wave function has the oscillation 
form 

Nn(x) = cos 2 <t>N x {x) + sin 2 <j)N 2 {x). (10) 



Perhaps more surprising, but unambiguously 
true from Eqs.(6) and (10), is the notion that the 
decay Eq.(3) of a muon with fixed four momen- 
tum throws the other particle 0, || (here the 
W^f f ) into an oscillation superposition of ampli- 
tudes |, | 

W+ f (x) = cos 2 <t>W+{x) + sin 2 (f>W+(x). (11) 

In Eq.(ll), W+ a cx (* 7a (l - 75)^) for j = 1,2. 

The four momentum kinematics may now be 
written as follows: Let be the initial four mo- 
mentum of the muon, and pj for j = 1,2 the four 
momentum of the possible neutrino mass eigen- 
states. Then 



Since total four momentum is conserved at each 
vertex of a Feynman diagram, even for internal 
virtual particles, one then associates two possible 
four momenta for the 

P eff w,j = (P^-pj), (i = 1,2) (13) 

whose effective (virtual off shell) mass M e ffw, 

(M effw c) 2 = -P* ffwA = -Peffwv (14) 

is much smaller than the on shell mass of the W + ; 
i.e. M eff « M w . From Eqs.(12), (13) and (14) 
it follows that 



-K-c) 2 , (j = l,2). (12) 



2(p 2 - Pl )-P f2 = {m 2 ~m 2 )c 2 



(15) 



In the rest frame of the muon, the two possible 
neutrino energies, ei and e 2 have a Bohr transi- 
tion frequency uj determined by Eq.(15); It is 



flUJ = 62 — £l — 



\ m l 2 
2M„ 



(16) 



The central result of this work now follows from 
Eq.(6) for the virtual wave function. This 

wave function contains a product of a neutrino 
wave function N^{x) and a muon wave function 
ifr(x). After absolute value squaring the W^j de- 
cay amplitude, one finds (from the neutrino wave 
function) the usual neutrino oscillation probabil- 
ity 



cos 2 6e~ l ' lt/h 



sin 



(l - sin 2 (2(f)) sin 2 (ut/ 2)), 



(17) 



and (from the muon wave function) the usual 
muon decay probability exp(— Tt). Thus, for the 
probability P(t) of the W^jj decay survival, we 
have the differential equation 

dP(t) 



dt 



ZTe~ rt (l 



(2(j})sin 2 (uit/2)).(lS) 



In Eq.(18), the constant Z is found by normaliz- 
ing the the total W^* survival probability; i.e. 



dP 
~dt 



dt = P(0) - P(oo) = P(0) = 1. (19) 



To test the above notion of observing oscil- 
lations in the muon neutrino P^, we analyzed 
the KARMEN anomaly data using our central 
Eq.(18). In the KARMEN experiment §, §, §, 
the electron neutrinos from /i + — > e + + za, + v e 
were detected at times {ti} (after the muons were 
stopped) within time bins of width Atj = 0.5/xs. 

One expects (using our theory) a mean number 
of detected electron neutrinos (from the muon de- 
cays) in each bin to be given by 



fJ,i{uJ,4>) = 



dti 



At, 



(20) 



where Af is chosen so that fa is the (total) 
number of observed v e events. If the events in 
each time bin obey Poisson statistics, then the 
probability distribution for a counting sequence 
of rii events in the i th bin is given by 



V\n\ uj, 1 



n 



n,\ 



(21) 



The likelihood function (for the theoretical pa- 
rameters uj and 4>) is then determined by experi- 
mental KARMEN data n, data via 



A(cj, <j>) = JCP[n data ; uj, 0] 



(22) 



where K, is an arbitrary constant. Under the 
best circumstances, the likelihood function ex- 
hibits a single maximum in the theoretical pa- 
rameter space (oj,4>). 

We employed this method of determining the 
theoretical parameters uj and </>, thus analyzing 
the KARMEN data from the viewpoint of the 




Figure 1. The likelihood as a function of the pa- 
rameters u> and <j). 



induced oscillating W^j decays. The likelihood 
function is exhibited in Fig. 1 . A single maximum 
likelihood peak was found. The derived values of 
uj and <p were determined to be 

u = (1.60±0.16)//xs, (23) 

and 

= 0.34 ±0.10, (24) 

within a 95% confidence interval obtained from 
the likelihood function. In terms of the neutrino 
mass squared differences, 

|™2-™i| = (0-22 ± 0.02) (eVyc 2 ) 2 , (25) 

also within a 95% confidence interval. 

In Fig. 2, we compare the theoretical real time 
oscillation in the v e counting rate (from the )i + 
decays) with the experimental oscillation for the 
v e counting rate in each time bin. The theoretical 
parameters (w, (f>) were chosen from the maximum 
of the likelihood function, The over all fit may be 
described by the statistical (x 2 /dof) ~ 1.4. 
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Figure 2. The theoretical time dependence of the 
v e counts along with the experimental points. 



In summary, the real time KARMEN anomaly 
for the electron neutrino counting rates (obtained 
from muon decay) has been analyzed. The data 
provides evidence in favor of a neutrino flavor ro- 
tation model in which the muon neutrino is in 
a superposition of two mass eigenstates. The 
neutrino mass splitting fits the data employing 
Eq.(25), and the flavor rotation angle fits the data 
employing Eq.(24). The theoretical fit to the ex- 
perimental KARMEN anomaly appears quite sat- 
isfactory. 
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